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Abstract: 

A grand-canonical system of interacting bosons is considered to study phase 
transitions of ultracold atoms in an optical lattice. The phase diagram is dis- 
cussed in terms of a matrix-like order parameter, representing a symmetric phase 
(Mott insulator) and a symmetry-broken phase (Bose-Einstein condensate). An 
exact solution was found by introducing N colors and taking the limit N — > oo. 

1 Introduction 

In this paper an interacting Bose gas in an optical lattice is studied. It is based 
on some previously published ideas for the hard-core Bose gas and its gener- 
alization to soluble models. The application to Bose-Einstein condensation, 
depletion of the condensate and the Mott insulating state in these systems is 
discussed. 

A real Bose gas, as has been studied in the case of trapped cold alkali atoms 
[Q, is always subject to interparticle interaction. In a good approximation this 
interaction is a hard core, characterized by a scattering length a. Depending 
on the type of atoms and external conditions (e.g. absence or presence of an 
external magnetic field), the scattering length can vary over a wide range. Typ- 
ical values of a observed in experiments are a ss 5nm in 87 Rb [2] and up to 
a « 500nm in a gas of 85 Rb atoms near the Feshbach resonance 'A . Negative 
scattering lengths are also known but will not be considered here. The effect of 
the hard-core interaction can be characterized by the interaction parameter 

j scattering length 

typical interparticle distance 

If the typical interparticle distance in a Bose gas is large (e.g., « 200nm in 
a very dilute gas that undergoes a Bose-Einstein condensation [I]), / is small 
(7 w 0.01), and the interaction can be treated as a weak perturbation. This is 
the regime where the Gross-Pitaevskii physics can be used 0J E] • 



If the Bose gas is transferred to an optical lattice the characteristic length 
scale is given by the distance of the potential minima, the lattice constant a;. 
Typical experimental values are a; w 250...500nm. Then the interaction is 
/« 0.01. ..1. 

^From the theoretical standpoint the simplest case is the non-interacting Bose 
gas (I — 0), where no qualitative effect of the lattice on the Bose-Einstein con- 
densation is expected. Indeed, Bose-Einstein condensation in a non-interacting 
Bose gas is not a phase transition, characterized by spontaneous symmetry 
breaking, but a macroscopic occupation of the groundstate. A grand-canonical 
ensemble of non-interacting bosons at inverse temperature /3, with quantum 
numbers I = 0, 1, ... and energy e; > is given by the partition function 

Z Q =Tre~^ = Y[Y J {Ce^ l f ni . (1) 

;>o n,>0 

The fugacity £ is chosen such that eo = and e; > for I > 0. It is assumed 
that the groundstate is unique. Then ( > 1 is not allowed here because this 
would lead to a divergent series in Eq. for I = 0. The average number of 
particles with quantum number I is 

This quantity decreases with increasing e/ and has its maximum at I = 0. The 
probablity of finding a particle in the groundstate is 



No =1 

N N 



The Bose-Einstein condensate (BEC) is the macroscopically occupied ground- 
state, i.e., lim7v-»oo N /N > 0. According to the result this can only happen 
if C goes to 1 when N goes to infinity. Using the density of states p(e) at energy 
e and keeping the ratio N/V (V: volume of the Bose gas) fixed when TV — > oo, 
the second term in Eq. becomes the integral 

v r^/\. 



N/3 Jo e x -l 

There exists a BEC if the integral is smaller than 1. The latter is the case if 
(i) the density of states p(e) cuts-off the pole of the integrand (this is usually 
the case for three-dimensional lattices) and if (ii) V/{N0) is sufficiently small. 
Thus the presence of a lattice does not prohibit the formation of a BEC in 
a non-interacting Bose gas at sufficiently low temperatures. More interesting 
physics is expected for J w 1 when the lattice constant is of the same order as 
the scattering length because the interaction can block the tunneling processes 
and lead to an insulating state without phase coherence and with a locally fixed 
number of bosons, the Mott insulator. 
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2 Random- Walk Representation 



In order to discuss the interaction between bosons in an optical lattice it is 
convenient to derive a random- walk representation for the lattice Bose gas. This 
is related to the idea that the average over particle fluctuations in an equilibrium 
system is equivalent to the average of random walks. For this purpose one can 
start from non-interacting bosons, then develop the weakly (Bose-Hubbard-like) 
and finally the strongly (hard-core) interacting gas. 

2.1 Non-interacting Bose Gas 

Evaluation of the geometric series in the partition function of Eq. gives 

1 



Z n = 



fir 



Assuming that (3 is an integer, the identity 

l-(f e -f*i = det(l - Cu»i) 

with the j3 x (3 matrix 

/° 
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(4) 



and the diagonal (3N x j3N matrix w 
write the partition function as 



Z 



det(l - Cw) 



diag(wi,W2, ■■■,wat) can be used to 

(5) 



The matrix elements in the last equation are defined with respect to x = (I, t) 
(t = 1, 2, I = 0,1, ...,Af). The diagonal matrix w can be replaced by any 
other matrix to describe a non-interacting Bose gas. For instance, tunneling 
between sites r, r' of the three-dimensional optical lattice gives with x — (r, t) 
the matrix elements 

{JSt^t+i/6 for r, r' nearest neighbors 
(1 - J)5 t ',t+i for r' = r 
otherwise 

This choice implies that a particle either stays in a potential well of the optical 
lattice at r with probability 1 — J (0 < J < 1) or tunnels from a well at r to 
one of the six neighboring wells with probability J/6. An expansion of Z a in 
powers of C 



(6) 



n>0 



(7) 
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has an interpretation in terms of random walks: the matrix elements w XjX > can 
be considered as (random) steps in positive time direction between the lattice 
sites at r and r'. Then the trace Tr(w° l ) consists of terms which start at a site 
r, go 1(3 steps and eventually return to r. This is a consequence of the fact that 
w has non-zero matrix elements only with t,t + 1. For I > 1 the walk is sent 
back to t = 1 after (3 steps due to the periodic boundary conditions defined by 
the matrix A realization of random walks is shown in Fig.l. The number 
I — 1 can be associated with the winding number of the corresponding walk. 



2.2 Bose Gas with Intermediate Interaction 

A weakly-interacting Bose Gas is described by the Bose- Hubbard model 6 
[3 This can be derived from the random- walk expansion of Eq. by 
adding an interaction between the random walks. Since the interaction is local 
and repulsive, the statistical weight of intersecting walks must be reduced in 
comparison with their non-intersecting counterparts. This can be implemented 
in the random-walk representation of Eq. {JJ by introducing a phase factor to 
the tunneling matrix element 

^x,x' * ^x/x'^- ^ = V x x '. 

and using for the partition function of the interacting system 

Zi = (detil-Cv)- 1 )^ 

with the Gaussian distribution of ip x with zero mean, and {^p x )ip = U. The 
expansion of Zi in powers of £ yields terms of the type 

x' x' 

For a singly occupied x (i.e. ^2 , l x i = 1) the averaging gives a rescaling of the 
fugacity £ — > (,e~~ u . In the case of a single well (i.e., J = in expression 
the partition function can be easily evaluated and gives with £ = e M (ji being 
the chemical potential) 

" e -pu(n-n/u)n^ 



n>0 



The summation over n > 0, here formally obtained from a geometric series, 
corresponds with the summation over particle numbers, as one can see by com- 
paring with Eq. IjTJl. For large (3 only the term with minimal value of (n — n/U)n 
contributes significantly. This is a (discrete) saddle point approximation: 



e -PU{n-^/U)n ^ maxe -PU(n-n/U)n 

Z-~i ^ n>0 

n>0 
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The value of n that gives the maximum is also the average number of particles. 
To evaluate this quantity one has to determine the groundstate energy 

E (jjl/U) =minU(n- n/U)n 

for a given value of fi/U. The energy E(/j,/U,n) = U(n — fi/U)n as a function 
of /i/U is linear. This implies for the groundstate energy a piecewise linear 
function with slope — n for 2n ~ 1 < fi/U < 2n+ 1 (cf. Fig. 2). If one includes 
tunneling processes, e.g. using an expansion in small J, the phase diagram 
would give the lobe structure, discussed for the Bose-Hubbard model |SJ[71|S]. 



2.3 Hard-core Bose Gas (/ = 1) 

In the case of hard-core interaction it is not just that random walks have to pay 
a penalty when they intersect but intersections are prohibited completely. This 
is a constraint to the random walks which can be included explicitly by using 
the functional integral representation. Starting from Zq 



1 



det(l - Cw) 



J exp ( - £ - :,,'!-;.) Yl d$ x d%/>ir, 



the complex field $ x is replaced by a nilpotent field: & x — > rj^, $* — > 77^, where 
(viY = for Z > 1. With the integral over the nilpotent field 

TT TT 77^ = / 1 if^ / = t,j s = l,2 j 
h c t / t \- X 1 otherwise 

where L is the full lattice spanned by x, the partition function of the hard-core 
Bose gas reads 



Zhc = exp [ V ^(1 + (w),^^, 

J HC 1 — ' 



The expansion of Zhc i n powers of C produces the non-intersecting random 
walks thanks to the nilpotent condition of the 77's. 

Physical interesting quantities like the density of particles n are expressed 
directly in terms of Zhc'- 

C d log Z HC ^ 



(3Af dC, 



The groundstate of a Bose gas can exhibit off-diagonal long-range order that 
is characterized by a non- vanishing condensate density no, obtained from a 
correlation function that is related to Zhc '■ 



I™ y— -is J2 I r & r fir,t) ex P [53^(1 + (w)x,x>vl 
IrHoo Zhc p ^r 1 Jhc 



t JHC x,x> 



= n Q > 0. 
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3 Fermion Representation of Hard-Core Bosons 

The nilpotent variables rf x can be replaced by products of Grassmann variables 
as rfj. — > ( — l^ipx^x- This gives with the usual Grassmann integration [5] the 
partition function |14| 



with 

x x,x' 

The advantage of this representation is that a quadratic form of Grassmannians 
in the exponent can be integrated out, leading to a determinant. Of course, 
in the case under consideration Sf erm is not a quadratic form. However, the 
non-quadratic part can be expressed by a coupling of a quadratic form to an 
Ising spin S x ,x' = ± 1; 



exp (C Wx ^' ^xV>x' 



= 2-^ J2 ^(E^^^+^J' 

{S x-X ,=±l} 

where Af is the number of lattice sites. Now the Grassmann integration can be 
performed and gives 

Z HC = 2-^ J £ det(l + U) 2 

{S x , x ,=±l} 

with the matrix elements 



U Xl x' — \ C l ^x^x' S x ,x' • 

The density of bosons n is the expectation value 

1 £ { s.., } det(l + U) a Tr((l + U)- 1 U) 

77 = ■ 

W E {S ^ } det(l + U) 2 

^<Tr((l + U)-U) )s , 
with respect to the distribution of Ising spins 

dct(l + U) 2 



P(U) 



This representation is a good starting point for a Monte Carlo simulation of the 
hard-core Bose gas 112) . 
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In one dimension the partition function can be evaluated exactly, since pe- 
riodic boundary conditions are satisfied for each random walk individually (i.e. 
the winding number is zero) due to the hard-core condition. Then the determi- 
nant is just a sum over all products along possible random walks with periodic 
boundary conditions in t: 

det(l + U) = 1+ X) II 

ir=jtid r 

with 

9r,iT — II ^c,7r(a:) • 
t 

Each g r ^ contains a product of Ising spins. Now one takes the square of the 
determinant and the sum with respect to the configurations of the Ising spins. 
This gives 

e (i+eii^) 2 =2^+e e n^^^+Ena 

{S xx ,=±l} ir^id r njHd{S XiX ,=±l} r ix^id r 

since does not depend on the Ising spins: 

9r,n = Y[ U x,n(x) = C P Ylw xMx) . 

t t 

Therefore, the partition function can also be written as a determinat 

Z HC =det(l + Cw). 

For the occupation of the quantum number I follows 

_ Cfe^ 1 f 1 for Ce e ' > 1 

1 ~ 1 + C,&e^ ~ \ for Ce £i < 1 

for [3 ~ oo. Thus a quantum state is either empty or singly occupied, a conse- 
quence of the hard-core interaction. This is also the property of a non-interacting 
fermion gas. It reflects the equivalence of the hard-core Bose gas and the free 
fermion gas in one dimension |13j . 

4 Generalization to N States and N Colors 

Generalizations of the hard-core Bose gas have been proposed that lead to ex- 
actly soluble limits: the N state bosons (NSB) ^2] and the N color bosons 
(NCB) ^H]. In the first case a boson can occupy one of N different states at 
each lattice site, where these states have the same energy. The hard-core inter- 
action applies only to bosons in the same state. Therefore, bosons can avoid 
the interaction during the tunneling process by choosing an unoccupied state. 
If N is large the NSB are effectively weakly interacting. Nevertheless, it turns 
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out that the interaction plays an important role even in the limit N — ► oo |14| . 
The other generalization uses colored bosons: the hard-core interaction applies 
only to bosons with the same color. In this situation the interaction is effec- 
tively strong for all values of N. In particular, it will be shown that even in the 
limit N — ► oo it can destroy the Bose-Einstein condensate and create a Mott 
insulator, an effect not possible in the limit N — > oo of the NSB. 



4.1 Effective Field Theory of NSB 

The extension of Sf erm to the NSB is given by 



N N 



snsb = E E(^ ,a ^' Q + ^l' a ) + ^ E w ^' E ^l-^/tif- 

x a — 1 x,x* a,/3=l 

The Hubbard-Stratonovich transformation 

E^' a €' Q ^^,x, 

a 

leads to the partition function 

Znsb = J e- NS '»s*V[<p, X ] 

with the action 

S' NSB = (cp, (1 + w)-V) + (x, X*) - E 10 S[C _1 + + + ixl)}- (9) 

X 

S' NSB does not depend on N and has a U(l) symmetry: it is invariant under 
the transformation ip x — > e lr, <p x , Xx —> e tv Xx, with a global phase rj. 



4.2 Effective Field Theory of NCB 

In the case of NCB the action Sf erm is extended to 

x a—1 x,x' a,/3=l 

The Hubbard-Stratonovich transformation 

a 

leads to the partition function 

Z NCB = f exp{ - TVE +2iVlogdet(l + u)}p[u]. (10) 
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With the choice of w given in Eq. © the action S can also be expressed, after 
a decomposition of the matrix field u = uy + uj^ (uy is diagonal in space), as 

S'ncb = ^rTj) Tr ( U ll U f) + ^ Tr ( u ^ u ^) - 21ogdet(l + U|, + uj.) 
This is invariant under the orthogonal transformation 

u -> OuO T . 



5 TV — > oo Limits 

The new representations of the partition functions depend on N only through 
a prefactor in front of the TV-independent action. Therefore, the functional 
integration can be performed for large values of N by the saddle-point (SP) 
approximation. That means we have to solve the extremal condition of the 
action S: SS = 0. Systematic corrections to the SP solution are obtained in a 
1/N expansion. 



5.1 Results for the NSB 

There are two SP solutions, one is the trivial solution tp = x = which is valid 
for C < 1 an d a non-trivial solution with \\\ = |<p|/2 = yl — 1/C which is valid 
for £ > 1. This implies for total density of bosons at zero temperature 

,fs JO for C < 1 

~ \ 1 _ £-1 for C > 1 ' 

Moreover, all bosons are in the condensate in the N — > oo limit. However, the 
1/N corrections indicate a depletion of the condensate, as shown in Fig. 3. This 
phenomenon reflects a partial destruction of the condensate due to interaction. 
It can be understood as a precursor for a new groundstate in the interacting 
Bose gas, visible in the NCB as the Mott insulator (s. next section). The 
critical point £ = 1 of the Bose-Einstein condensation does not depend on the 
tunneling rate J and agrees with J = for all N. It should be noticed that the 
limit N — > oo of the NSB leads to the Gross-Pitaevskii equation JU], at least 
in a close vicinity of the critical point (dilute regime). 

Quasiparticles are available from 1/N corrections. In the condensate (i.e. 
C > 1) they have a dispersion on large scales with respect to the three-dimensional 
wavevector k as 

e(k) = y j(C ~ 1) fc2 + (j/ 6 )2 fc 4 

Due to the spontaneously broken U(l) symmetry of the NSB there is a Goldstone 
mode (f> for £ > 1 which satisfies the equation 



[c 2 k 2 -uj 2 ^)4>(k,uj) = 0, 



where c = y J(£ — 1) /3 is the sound velocity. 
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5.2 Results for the NCB 



For the SP equation of the action Sncb one can distinguish symmetric and 
symmetry-breaking solutions. Examples of the former are the trivial solution 
Up = Ui = (empty phase) and Uy ^ 0, = with identical matrix elements 

_ / My for r' = r, t' = t + 1 

X X ' — 1 * 

L otherwise 



This ansatz leads to the mean-field action 

Sncb = M(i 



2 

llog(l + <) 



C(l-J) /3^^"\ 
which is minimized by the solution of the SP equation 

For zero temperature the value of u\\ jumps from to 1 at £ = 1/(1 — J). The 
density of bosons can be calculated with Eq. (JSJ) from the mean-field action 
and gives n — 0(C(1 — J) — !)• This is the behavior of the system without 
tunneling. The jump should be intercepted by an intermediate phase, most 
likely the Bose-Einstein condensate. The off-diagonal long-range order of the 
NCB in the large-TV limit is related to a symmetry-breaking SP solution with 
7^ 0. This represents a Bose-Einstein condensate with a Goldstone mode. It 
is anticipated that this solution agrees qualitatively with the result of the NSB, 
discussed in the previous section. The resulting phase diagram, obtained from 
the N — > oo limits of the NSB and the NCB, is shown in Fig. 4. 



6 Conclusions 

The hard-core Bose gas at zero temperature has been studied in terms of the 
soluble N — -> oo limits of NSB and NCB. NCB describe three different phases 
of the hard-core Bose gas, characterized by a matrix-like order parameter u = 
Uy + u^. There is an empty phase (uy = = 0), a Bose-Einstein condensate 
(uy = 0, uj_ 7^ 0), and a Mott insulator (un ^ 0, = 0). In contrast to the 
NCB the complex scalar order parameter of the NSB describes only an empty 
phase and a Bose-Einstein condensate but not the Mott insulator. Thus, the 
NSB has only one symmetric SP solution (empty phase), whereas the NCB has 
at least two (empty phase and Mott insulator). However, the depletion of the 
condensate in o(l/N) of the NSB indicates the partial destruction of the latter. 
The main result of this study is the effective action of the NCB of the matrix 
field u 

S'ncb = m3j) Tr ( u ll u f) + lJ Tr ( u ^ - 2io Sdct(l + u,| + u x ), 
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where its minima represent different phases of the strongly interacting Bose gas. 
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Figure 1 : Contribution of the random- walk expansion to the partition function. 
Sites a and b are connected by the periodic boundary conditions, a consequence 
of the fact that the partition function is a trace (cf. text). The hard-core 
interaction does not allow intersecting random walks. 



12 




3 4 
CHEMICAL POTENTIAL 



Figure 2: Groundstate energy and step- like density of particles of an interacting 
Bosc gas in a single well as a function of n/U. 
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Figure 3: l/N corrections of N state bosons: enhancement of the total density 
n (upper three curves) and depletion of the condensate density no (lower three 
curves) [TS] . 
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Figure 4: T = 0-phasc diagram obtained from N — > oo limits. The phase 
boundary to the Mott insulator was calculated with N color bosons, the phase 
boundary to the empty phase with N state bosons. 
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